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ABSTRACT 

We use a volume-limited spectroscopic sample of isolated galaxies in the Sloan Digital Sky Survey (SDSS) 
to investigate the frequency and radial distribution of luminous (M r < —18.3) satellites like the Large Mag- 
ellanic Cloud (LMC) around ~ Milky Way analogs and compare our results object-by-object to ACDM 
predictions based on abundance matching in simulations. We show that 12% of Milky Way-like galaxies host 
an LMC-like satellite within 75 kpc (projected), and 42% within 250 kpc (projected). This implies ~ 10% 
have a satellite within the distance of the LMC, and ~ 40% of galaxies host a bright satellite within the 
virialized extent of their dark matter halos. Remarkably, the simulation reproduces the observed frequency, 
radial dependence, velocity distribution, and luminosity function of observed secondaries exceptionally well, 
suggesting that ACDM provides an accurate reproduction of the observed Universe to galaxies as faint as 
L ~ 1O 9 L on ~ 50 kpc scales. When stacked, the observed projected pairwise velocity dispersion of these 
satellites is a ~ 160 kms -1 , in agreement with abundance-matching expectations for their host halo masses. 
Finally, bright satellites around primaries are significantly redder than typical galaxies in their luminosity 
range, indicating that environmental quenching is operating within galaxy-size dark matter halos that typically 
contain only a single bright satellite. This redness trend is in stark contrast to the Milky Way's LMC, which is 
unusually blue even for a field galaxy. We suggest that the LMC's discrepant color might be further evidence 
that it is undergoing a triggered star-formation event upon first infall. 

Subject headings: galaxies: Local Group — galaxies: Large Magellanic Cloud — galaxies: evolution — 
galaxies: dwarf — galaxies: groups: general — cosmology: dark matter — cosmology: 
observations 



1. INTRODUCTION 

In the ACDM concordance cosmology, dark matter sub- 
structure is ubiquitous due to the continuation of the power 
spectrum to very small scales. As a result, substructure pro- 
vides valuable constraints both on the nature of dark matter 
and its effects on galaxy formation. The missing satellites 
problem — the observation that ~ galaxies like the Milky 
Way (MW) have too few satellites for the predicted amount 



of dar k matter substructure (Klypin et al. 1999; Moore et al. 
1 1999)1 — has been discussed in depth for taint satellites. Pro 



posed solutions incl ude steep slopes or truncations in galaxy 
formation (e.g.|Bullock et al . |2000| | S trigari et al.|2008||Bovill| 
|&Ricotti|20 09; Kop osov et al.l2009||Kjavtsov|2010| and ref- 
erences there in), as well as a si gnificant component of obser- 
vational bia s UTollerud et al.|2008||WaTsh et al.|2009||Bullock| 
|et al.|2Q10] >. 

In contrast, less is known about more luminous satellites of 
~ hosts, e.g. satellites like the Large Magellanic Cloud 
(LMC) or Small Magellanic Cloud (SMC). These luminosity 
ranges are instructive to compare to ACDM, as the substruc- 
ture is significant for ~ hosts at these scales, while the 
satellites are still bright enough to be detected at cosmological 
distances. Thus, unlike faint Galactic satellites, these brighter 
satellites can be studied outsid e of the Local Group (LG) 
with statisticall y large samples ( Woods et al. 2006| |Woods| 
& Geller|2007| l. Studying the small-scale structure (the halo 
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substructure at cosmological distances) of isolated galaxy ha- 
los therefore provides a valuable test for ACDM. 

In ACDM, L MC-sized satellites may be rare around MW- 
like hosts (e.g., IBenson et al. 12002 [|S trigari et al.12007) |Ko-| 
posov et al. 2009| l. |Boyl an-Kolchin et al. (2010~b| consid- 



ered directly how common MW/LMC and MW/LMC/SMC 
systems are in a very high-resolution ACDM simulation (the 
same discussed in Sj3] below), showing they are potentially 
quite uncommon (3 — 25% ), depending on the assumed MW 
mass. |Busha et~aL (2010b} performed a similar analysis on a 
different simulation, with similar result^] 

Additionally, there are observatio nal hints that LMCs are 
unusual around MW-like hosts (|James et al. i 2008 ; James & 
[lyor^[20T0t [Liu et al |[20T0l |van den Bergh|]201o| . Despite 
this possibility, models are sometimes tuned to produce the 
correct number of bright satellites like the MCs, or are ad- 
vanced because they naturally do so (e.g. ILibeskind et al. 



2007| |Koposov et^|2009t |Kravts ov 2010). Thus, an un- 
derstanding of the abundance and distribution of luminous 
satellites around isolated hosts constrains models of galaxy 
formation, and helps place the Local Group itself in a clearer 
cosmological context. 

In this study, we present a set of spectroscopically selected 
galaxy pairs in the Sloan Digital Sky Survey (SDSS) com- 
posed of an isolated ~ primary and an ~ £lmc secondary. 
We statistically connect these systems to ACDM dark matter 
halos selected identically from simulations. While we would 
ideally include SMC-like satellites, the SDSS spectroscopic 
survey is not deep enough to provide a statistically significant 

4 While these results, and ours, may show differences at a factor of two 
level, we show below that this is primarily due to the exact choice of defini- 
tion for "MW-like" or "LMC-like" and thus all are in accord when compared 
appropriately. 
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sample of SMC analogs. We thus probe the clustering prop- 
erties of the lowest-luminosity satellites possible with SDSS 
spectroscopy, and thereby provide a robust test of ACDM for 
the smallest subhalos and scales that can be probed with a 
statistical significant sample. Given the similarity of our sys- 
tems to the MW/LMC group, this approach also allows us to 
address how common such systems are in the local universe. 

Throughout this paper we use the term "primary" and "sec- 
ondary" to indicate the larger and smaller galaxies in our pro- 
jected pair samples, for which the pair may not actually be 
physically associated. In contrast, we use the terms "host" 
and "satellite" to indicate pairs for which the smaller galaxy 
is definitely within the larger's halo (and hence only applies 
for halos in the simulation, where full 3D information is avail- 
able). 

This paper is organized as follows: In *]2jwe describe how 
we construct our observational sample, inH5]we describe the 
simulation we use to interpret this sample, in ffflwe analyze 
the sample and compare to the expectations of ACDM, in £j5] 
we discuss the MW/LMC system in the context of this sample, 
in f|6]we compare our results to recent related works, and in 
f|7]we review our primary results. 

2. OBSERVATIONAL DATA 

We begin by describing the SDSS data sample chosen to 
identify potential host/nearest satellite pairs (i.e., primaries 
and secondaries) in a volume-limited sample . We use the 
NYU V alue-Added Galactic C atalog (VAGC |Blanton et al. | 
|2005b| > based on SDSS DR7 ( jAbazajian et al.||2009|l with" 
the im proved photometric calibrations (Padmanabhan et al. 
2008). From this cata log, we examine the main galaxy sam- 
ple, using K-corrected Petrosian ( 1976 1 magnitudes with z = 



band shift. We use the r-band, which typically has the high- 
est signal-to-noise ratio of the SDSS bands. Where relevant, 
we assume WMAP7 (Komats u et al.|2010| l cosmological pa- 
rameters with h = Ho/100 = 0.704. 

We apply the following set of selection criteria in order to 
construct a clean sample of primaries and secondaries. As 
summarized in Table [T] we also construct samples using sub- 
sets of these criteria in order to test our method and to estimate 
completeness in our spectroscopic sample. 

1- A/r,pri - 51og(7i) < -20.0 (i.e., M r<pli < -20.8) for 
primaries, 

2. -20.0 < M r , sec - 51og(/i) < -17.5 (-20.8 < 



M r 



< —18.3) for secondaries, 



3. z pr i < 0.0342, which ensures the sample is complete 
based on SDSS spectroscopic survey completeness lim- 
its for M r - 5 log(/i) < -17.5 (M r < -18.3). 

4. The primary must have a measured redshift. 

5. The primary must be isolated: at most one other galaxy 



proj 



< 



in the primary magnitude range within 250 < d 
700 kpc/h (355 < d proj < 995 kpc) and with redshifts 
such that Av pr [ < 1000 kms" 1 . Ther e must be none 
within 250 kpc/h (355 kpc). (see p3|l 



6. The secondary must be within 250 kpc/h (355 kpc) 
projected and be the nearest (projected) galaxy in the 
secondary magnitude range. 

7. The secondary must have a measured redshift. 
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FIG. 1. — Examples of SDSS primary/secondary pairs in the clean sample 
(upper) and false pairs (lower). Secondaries identified by our criteria (see 
text) are marked with red circles (upper panels) or magenta triangles (lower 
panels). The upper three are all in the clean sample (have redshifts close to 
the primary) and span a range of projected separations. For the lower three 
images, blue circles are SDSS pipeline photometric objects, clearly showing 
the identification of HII regions as photometric objects. For these same lower 
three, the secondaries are clearly HII regions in the primary (or satellites that 
are indistinguishable from HII regions). We visually identify and remove all 
pairs of this kind from our sample. 



8. Av = c\ 



< 500 km s 



We now consider these criteria and their motivation in more 
detail. We only consider galaxies with measured redshifts as 
potential primaries. With these redshifts, we select a volume- 
limited sample of galaxies, with the maximum redshift set 
such that all primaries and secondaries are within the com- 
pleteness limit of the SDSS spectroscopic survey (r < 17.77). 
From these objects, we identify all galaxies more luminous 
than our primary/secondary boundary (M r < —20.8) as pri- 
maries. This sample has a median redshift of z mc d = 0.028, 
near the upper limit due to the larger volume at higher red- 
shifts. 

Without further selection, many of these primaries are in 
groups or clusters. In such dense environments, secondaries 
with even small projected separations from the primary are 
often not true satellites (see Sj3}. We therefore apply addi- 
tional criteria to the primaries to ensure we are sel ecting rela- 
tively isolated primaries. Specifically, following Bart on et al.| 
( 2007 ), for each primary we identify all other primaries within 
a projected distance of 1 Mpc for which the redshift difference 
between the two satisfies Aw pr j = c\z\ — Z2\ < 1000 kms -1 . 
An isolated primary is then defined as having at most one pri- 
mary within this distance, and if such a companion is present, 
it cannot be wit hin 3 55 kpc. As estimated using cosmological 
simulations in £ 3.3 this criterion results in ~ 90% of the sec- 
ondaries being true satellites, or a contamination of ~ 10%, 
while still maintaining a large sample. 

With a sample of isolated primaries in hand, we then iden- 
tify the nearest galaxy fainter than this boundary, but more 
luminous than the limit for secondaries (—20.8 < M r < 
— 18.3), and label this the secondary if it is within 355 kpc 
projected. For this step we assume that the redshift of the 
secondary is that of the host to determine absolute magni- 
tudes and physical distances. This defines our "Full" sam- 
ple of primary/secondary galaxies. To clean this sample of 
foreground/background galaxies that only appear associated 
with the primary in projection, we filter this sample such that 
the redshift of the secondary be measured (the "any z" sam- 
ple) and also require the secondary's redshift to be within 
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TABLE 1 

Summary of catalog samples. 



Sample Name 1 


Mag/Vol Limit 2 


Pri/Host Present 3 


Sec/Sat Present 4 


Sec Has Redshift 5 


At; Cut 6 


Obs or Sim 7 


N 8 


Primaries 


y 


y 


y or n 


N/A 


N/A 


Obs 


1075 


Hosts 


y 


y 


y or n 


N/A 


N/A 


Sim 


373477 


Clean 


y 


y 


y 


y 


y 


Obs 


467 


Simulated Clean 


y 


y 


y 


N/A 


y 


Sim 


149906 


Any z 


y 


y 


y 


y 


y or n 


Obs 


584 


No 2 


y 


y 


y 


n 


N/A 


Obs 


189 


Full 


y 


y 


y 


y or n 


N/A 


Obs 


773 


Control 


y 


n 


y 


y 


N/A 


Obs 


24298 



1 Name of the sample. 

2 Whether or not magnitude and volume limits are applied (Criteria |l|4) . 

3 Whether or not a primary (host in simulations) is present. 

4 Whether or not a secondary (satellite in simulations) is present (Criterion 6| 
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f, 



-2ll 



-20 



Whether or not a secondary must have a redshift (Criterion^. 

Whether or not a secondary (satellite in simulations) must be within 500 km s - 1 of its primary/host (Criterion[8j. 

7 Specifies if the data are from SDSS/VAGC galaxy catalog or the MS-II mock galaxy/halo catalogs. 

8 Number of objects in the sample. 

500 km s -1 of the primary. This defines our main sample of 
primary/secondary pairs (hereafter referred to as the "clean" 
sample), three of which we show in the upper panel of Figure 
[T] The remaining samples are determined by not including 
various steps of this procedure, and are specified in Table [T] 

One final filtering step is necessary due to imperfections in 
the SDSS photometric pipeline. For nearby bright galaxies (a 
significant portion of our primary sample), individual HII re- 
gions in a star forming disk are sometimes classified as sepa- 
rate objects in the SDSS pipeline, as shown in the lower panel 
of Figure [T] Because these HII regions are extended sources 
and have redshifts matching the host, they are included in all 
secondary samples, including the clean sample. Fortunately, 
this effect is only significant out to separations of ~ 10 kpc, 
so there are a small enough number of possible such pairs 
(~ 50) that they can be removed by visual examination. All 
samples below have been cleaned of these HII regions. This 
process is not ideal, as it may remove faint satellites that truly 
are present, but are superimposed over a disk and thus appear 
indistinguishable from an HII region. Alternatively, some 
secondaries may be hidden from view by being behind the 
primary. Thus, removing apparent HII regions may result in 
under-counting satellites at small projected separations. For- 
tunately, this is only of consequence for (f pro j S 10 kpc, rep- 
resenting only ~ 0.3% of a typical primary's dark matter halo 
volume. Hence, these under-counting effects are not signifi- 
cant relative to Poisson errors in the discussion below. 

The clean selection criteria define a sample of galaxy pairs 
that are isolated from other luminous galaxies, and com- 
posed of a ~ primary with a nearest satellite with L ~ 
Llmc- We show the g — r, M r color-magnitude diagram 
(CMD) as Figure [2] with the clean sample represented by 
the red (primary) and white (secondary) points plotted over 
all isolated primaries (above the black line) and all galax- 
ies with secondary-like magnitudes (below) in the volume- 
limited sample. We note here that because the luminosity 
function climbs sharply for fainter galaxies, the typical lu- 
minosity of the clean sample is close to the faint cutoff for 
these samples: for the isolated primaries, the median a bsolute 
magnitude is M r = —21.3 (very close to M* from Blanton 
|et al.|2 003 ), while for the secondaries, it is M r = —19.3. It is 
important to also note that the secondary sample is not com- 
plete because the SDSS spectroscopic survey is not complete. 
The two main sources of systematic incompleteness are due 
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FIG. 2. — Color-magnitude diagram (g — r and M r ) for the SDSS data set. 
Over-plotted circles are the clean sample: primaries with secondaries are col- 
ored red (above the horizontal black line), and their secondaries (below) are 
white. The underlying hex-binned plot is a Hess-like diagram of the number 
of galaxies in each bin. Above the separation line, it shows the counts of all 
objects in the primary sample (isolated and inside the volume limit). Below, 
the bins are counts of objects in the Control sample (secondaries in the vol- 
ume limit). Thus, the discontinuity in the bins at the line is due to application 
of the isolation criterion above and not below. 



to two fiber collision effects. The first is collision between 
the primary and secondary, but this does not significantly bias 
our sample because the typical primary-secondary separation 
is larger than the on-sky fiber diameter for most of our objects. 
The second effect concerns the case where a potential satellite 
is near a bright background cluster, and hence is not available 
to be targeted for a spectrum in the SDSS sample. This effect 
is discussed further in Sj4] 
The selection criteria discussed above define a sample 
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roughly analogous to the MW/LMC system^] While there 
are 467 MW/LMC pairs in the clean sample, we find only 7 
objects if we instead use a volume limit appropriate for SMC- 
luminosity objects (~ 2 mag fainter than the LMC). As a re- 
sult, we only consider MW/LMC-like pairs for the remainder 
of this work. We defer further discussion specific to the LG 
to§] 

3. CATALOGS FROM THE MILLENNIUM II SIMULATION 

While the observational sample described in the previous 
sections is designed to produce a sample of isolated L* hosts 
with LMC-like satellites, the purity and completeness are not 
a priori known. Hence, we compare a high-resolution cos- 
mological simulation to our SDSS sample and search for 
host/satellite pairs, "observing" the simulation in the same 
way we observe the SDSS data ( |Barton et al. 2007). We pri- 
marily use this simulated catalog to make direct comparison 
with the SDSS sample, thus factoring out as many selection 
effects as possible in testing the theory. Moreover, these sim- 
ulated observations can be used to validate our observational 
sample selection by comparing the perfect position/velocity 
information in the simulation to our observational selection 
criteria. This step is crucial for informing the observations, as 
it would otherwise be nearly impossible to interpret the obser- 
vations due to the influence of projection effects and peculiar 
velocities. 

We use the high-resolution Millen nium II simulation (MS- 
II) of Boylan-Kolch in et al.| ( |2009| l as our source of dark 
matter halos with the assumption that all massive dark mat- 
ter halos host galaxies. The simulation resolves a volume of 
137 3 Mpc 3 , a nearly identical volume to that probed by our 
SDSS sample (z < 0.0342), evolved from z = 127 to 0. The 
background cosmology is: 



Kj-avts oVeral.|2004||Conroy & Wechsl er|2009{ |Moster etaL 



2010| |Guo et al.||20l0| l. The z; max values from the Simula- 
tion are computed for all subhalos at z = as either: (1) the 



fitot =1.0, On =0.25, fl b = 0.045, n A 
h = 0.73, er 8 = 0.9, n s = 1 , 



0.75, 



(1) 



where h is the reduced Hubble constant at redshift zero, erg 
is the rms amplitude of linear mass fluctuations in 8 h^ 1 Mpc 
spheres at z = 0, and n s is the spectral index of the primordial 
power spectrum. The force resolution (Plummer-equivalent 
softening) in the simulation is 1.4 kpc (comoving), and the 
mass resolution is 9.43 x 10 6 M Q . Dark matter structures 
are ide ntified by first usin g a Friends-of-Friends (FOF) algo- 
rithm (Dav is et al.||1985] > to ide ntify distinct halos, then ap- 
plying tlTe^uW™rJaigorithm ( JSpringel etal.|2001) to these 
halos to find substructure. Subhalos are also linked across 
simulation outputs in merger trees. The combination of these 
analyses enables us to reliably track redshift zero properties 
and histories of all subhalos with maximum circular velocities 
!)„,„ in excess of 30kms . For further details of the sim- 



ulation, halo catalogs, and merger trees, see Boylan-Kolchin 
|etaL|d^009> . 

3.1. Abundance Matching 

To connect our SDSS galaxies to MS-II halos, we require 
a mapping from galaxy luminosity to halo maximum circu- 
lar velocity v max . Here we use the well-established sub- 
halo abundance matching technique to do this mapping (e.g. 

5 While the presence of M3 1 may seem to imply the MW is not "isolated", 
the LG is in fact most likely composed of two isolated halos following the 
definition of isolation adopted here. Thus, if viewed from far away with 
typical redshift errors, the LG would be statistically impossible to distinguish 
from two more widely separated galaxies for most orientations. 



current value of w n 



if the subhalo is dominant in its FOF 



group, or; (2) the value of v max the subhalo had at the time its 
mass was at a maximum, if the subhalo is non-dominant. Lu- 
minosities are then matched to maximum circular velocities 
by setting n(> L) = n(> v max ). 

Matching the halo catalogs from the MS -II to the galaxy lu- 
minosity function of Blan ton et al.| ( |2003] >, we find our lower 
limit for selection of host halos is v max > 166.5 km s _1 from 
-Wr,host < —20.8, and the constraint for satellite (sub)halos is 
166.5 > u max > 94.8 krns" 1 , from M r<hoat < -18.3. Addi- 
tionally, given the median observed M r values of —21.3 and 
— 19.3 for the primary and secondary samples (see Sj2]), with 
abundance matching we can estimate the median (sub)halo 
w m „ for our galaxies. For primaries and secondaries, this is 
211 km/s and 112 km/s, respectively. These are both within 1 
km/s of the median v max in the Simulated Clean sample, im- 
plying our abundance matching results are consistent between 
our subsamples. 

We note that by using the parameters for the selection crite- 
ria directly from the mock catalog, we are ignoring the effects 
of observational errors in the luminosity function. To quan- 
tify the effects of these errors, we perform 1000 Monte Carlo 
simulations of the luminosity funct ion by sa mpling each bin 
from the luminosity function of Blanton et al. (2005a) assum- 
ing Gaussian errors. For each of these, we compute the im- 
plied v max using the abundance-matching scheme described 
above. For both the primary and secondary limits, the scatter 
about the median v maK is ps 0.3%. 

Beyond the observational uncertainty in the luminosity 
function, there is an intrinsic galaxy-to-galaxy scatter in abun- 
dance matching due to stochastic effect s in galaxy form ation. 
This scatter has been quantified in Behroozi et al. (2010]), but 
their analysis is for matching halo mass to stellar mass, where 
the systematic errors in estimating stellar mass are the domi- 
nant effect. Because we use luminosity instead of stellar mass, 
observational systematics are less important, but the results 
are potentially color-dependent. Additionally, even random 
scatter in abundance matching will result in a systematic shift 
of our v max limits inferred from our M r limits. This is be- 
cause the mass and luminosity functions are both rising to the 
lower mass/faint end, so scatter in abundance matching will 
shift more galaxies from below to above the limit than vice 
versa. This has the effect of systematically biasing the v max 
limit for a given M r . 

To estimate the magnitude of this bias in the v max limit for 
a given M r for our data set, we use the results of the cosmo- 
logical smoothed par ticle hydrodynamics (SPH) simulations 
of Si mha et al.| ( |2010[ ). They find that the la log-scatter of the 
ratio of the luminosity of a galaxy to the luminosity assigned 
by abundance matching is 0.15 dex. For each of the MS-II 
halos, we assign a median luminosity by abundance match- 
ing to the Blant on et al.| ( |2005aj ) luminosity function. We then 
draw 1000 samples from a lognormal distribution about this 
median with a = 0.15 dex. We use this ensemble of halo 
luminosities to produce 1000 simulated luminosity functions, 
and abundance match them to the halo w max function. Finally, 
for each abundance matching mapping, we determine the re- 
sulting f max limits for our primary and secondary cutoffs of 
M r < -20.8 and M r < -18.3. We find that the result- 
ing v max limits are indeed systematically biased low, but by 
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< 1%, which is much smaller than the scatter due to different 
lines-of-sight in MS-II (an estimate of cosmic variance). To 
produce scatter comparable to this cosmic variance effect, the 
scatter in luminosity must be a fact or of ~ 10, an implaus ibly 
large scatter based on the results of Behroozi et al. ( 20 1 0) > . 



l.Oi 



3.2. Mock Galaxy Catalogs 

With the delineation between hosts and satellites in the sim- 
ulation determined (hosts: u max > 166.5 kms~ ; satellites: 
Umax > 94.8 kms -1 ), we proceed to c reate mock observa - 
tions of the MS-II in a fashion similar to |Barton et al.| ( [2007| >. 
As described in f|2] the key elements of the observational sam- 
ple are a maximum distance for the volume limit (and as- 
sociated minimum luminosity), a lower luminosity limit for 
primaries, a lower luminosity limit for secondaries, isolation 
criteria for the primaries, and a cut on the difference in line- 
of-sight velocity between the primary and secondary. Table 
[T] specifies the exact criteria for the specific samples. Here 
we describe specifically how this sample is extracted from the 
MS-II catalogs. 

A random position within the box is chosen for the "ob- 
server", and a random rotation to the simulation volume is 
applied, placing the observer at the origin with the simula- 
tion volume in the positive (x, y, z) octant. Line-of-sight 
distances to all potential hosts in the simulation volume are 
computed, and for each of these, transverse distances to all 
other objects in the volume are calculated as well. All hosts 
within a transverse distance of 1 Mpc and with a line-of- 
sight velocity separation less than 1000 kms -1 of each host 
are then saved, as is the satellite with the smallest transverse 
distance. As in for observational data, we require (1) zero 
(at most one) additional hosts within 355 kpc (between 355 
and 1000 kpc) of each host; (2) z pri , z scc < 0.0342; (3) 

l^pri — z scc \ < 500 km s -1 . We replicate this procedure 250 
times in order to account for variations due to observer po- 
sitions. In later sections, we derive simulation quantities for 
each of these realizations separately, and use the median and 
68 % limits from this ensemble. 

In order to test the fraction of projected satellites that are 
true satellites - i.e., are physically associated with their host, 
and are not simply close in projection - we also compute the 
true distance to each identified satellite, and the true distance 
to the satellite that is closest in 3D (not in projection). These 
are often, but not always, the same (see the next sections). 

3.3. Isolation Criterion and Completeness 

In comparing our simulated galaxies to our SDSS sample 
we make use of a mock galaxy catalog that we select to expe- 
rience the same selection effects as the real data. In this sense, 
we do not need to worry about issues like completeness if our 
only aim is to test this particular ACDM model. Neverthe- 
less, we can use the simulations to estimate to what extent the 
selection criterion we impose actually select the galaxy halos 
we are attempting to select, and to determine the degree of 
completeness we obtain. 

For the SDSS sample, as described in fj2j we require a 
means of selecting primary galaxies to ensure they are iso- 
lated galaxies rather than in a group or galaxy cluster. To 



this end, we con sider an isolation criterion motivated by Bar- 
ton et al. (2007) such that the primary has no other primary 
within the radius in which we search for secondaries, and at 
most one out to a larger radius (criterion |5]from jpj. This cri- 
terion is intended to eliminate cluster galaxies andinteracting 
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FIG. 3. — Purity of secondary sample in MS-II. Each line shows the fraction 
of secondary that are actual satellites of the primary (subhalos in a FOF sense) 
as a function of projected separation. The solid (blue) line is for the simulated 
clean sample, green (dotted) is for the any z sample (no Av restriction). 
Dashed (red) is for a sample like the clean sample but without the isolation 
criteria applied to the primary, and dash-dotted (magenta) is any z without 
isolation. 

primaries, while keeping galaxies that are isolated or in very 
loose groups like the LG. 

Note that any isolation selection wil l necessitate a tr ade off 
between purity and completeness (e.g. |Liu et al.|2010} Figure 
1). Our primary concern is purity, and as illustrated below, 
our choice of sampl e selection does indeed perform well in 
that regard ( jBerrier et "aLl|2006| |Barton et aL"1|2007| |Berrier 



et al. 201 1 ). Purity is the goal because we aim to measure the 



fraction of primaries with a secondary. Therefore, the overall 
completeness is not important for the comparison of observa- 
tions to simulations if our incomplete selection is unbiased. 
Our sample does present one possible bi as, however: because 
ACDM features "assembly bias " ( |Gao et al |2005| |Wechsler 
et al.||2006 Croton et al.||2007"] ), sampling from isolated ha 



los in low density regions may result in satellite distributions 
offset from a more general sample. Fortunately, as we show 
later (S 4.2 1, any assembly bias that may be present because of 
our selection does not seem to affect the radial profiles of dark 
substructure in the simulations in a significant way. 

With a MS-II sample (the Simulated Clean sample) selected 
to match the SDSS sample, we can examine how effectively 
these isolation criteria function to provide a sample of real 
satellites. Figure [3] shows the purity: the fraction of pri- 
maries in MS-II for which the secondary is an actual satellite 
(a real subhalo of its host) as a function of projected primary- 
secondary distance. The different lines illustrate how the pu- 
rity changes for different sample selections. The solid (blue) 
line shows that our fiducial "clean" selection (isolation crite- 
ria plus Av condition) provides a sample that is quite pure, 
with 80 — 90% of identified secondaries being actual satel- 
lites. If the isolation criteria are dropped (dashed red) but the 
Av condition is kept, the purity plummets to 50 — 60%, sug- 
gesting that by selecting isolated primaries we are greatly en- 
hancing our ability to interpret the data. The dotted green line 
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shows that the Av cut is most important at large separations - 
even when we select primaries that are isolated by our crite- 
rion, without a Av condition the purity drops to ~ 65% at the 
250 kpc projected. As illustrated by the dot-dashed (magenta) 
line, less than half of the sample is pure in the absence of ei- 
ther criterion. In summary, our fiducial clean criteria provide 
a sample that is made up of ~ 85% true satellites within 250 
kpc. The remaining 15% are halos that are in the secondary 
mass range, but are isolated (field galaxies) and not actually 
subhalos of the primarjj^] 

A final question is the incompleteness due to the effect of 
projection in determining the closest bright satellite. Some 
fraction of the nearest neighbor galaxies in projection will not 
be the true closest satellites in 3D distance. While this is not 
important for comparing ACDM to observations, this does af- 
fect the observationally implied fraction of primaries with a 
true satellite. We can use the MS-II mock survey to estimate 
this effect in our simulated clean sample by simply comput- 
ing the fraction of primaries with a secondary within a given 
projected distance that have a true satellite (a subhalo of the 
appropriate mass) within the the primary's halo. Dividing the 
clean-sample purity (Figure[3]l within the same projected sep- 
aration by this completeness thus gives the actual fraction of 
true closest satellites implied by an observed fraction within 
that separation. The associated ratio is larger than 2 for any 
separation within 100 kpc, implying that the projected near- 
est neighbor is unlikely (less than 1/2 of the time) to be the 
true nearest neighbor at separations smaller than this. This is 
because the smallest line-of-sight separations tend to be dom- 
inated by chance orientation. Interestingly, the correction fac- 
tor is unity at a projected distance of ~ 250 kpc, implying 
that the observed fraction of nearest projected galaxies (in the 
clean sample) should be equal to the true fraction of nearest 
satellite galaxies (within the three-dimensional virial extent of 
the halo) at this distance. 

4. COMPARING THE OBSERVATIONS TO THE SIMULATIONS 

Having established the particulars of our observational sam- 
ple and identically selected comparisons in the MS-II simula- 
tion, we can consider whether the observed satellites are con- 
sistent with this (bare bones) ACDM model. We do this by 
comparing the observed galaxy and simulation halo samples 
through unambiguous measurements in both the observations 
and simulations: sky coordinates and line-of-sight velocities. 
We also compare luminosity functions for our observed sam- 
ple to those implied by abundance matching in the simulation. 
Additionally, we show this sample does not deviate strongly 
from the general halo population for ~ galaxies, and infer 
the 3D satellite distribution. 

4.1. Radial Profiles 

An important comparison between the observations and 
simulations is that of the projected spatial separation between 
the primaries and their secondaries (e.g ., |Chen et aL"1|2006[ 
Klypi n & Prada||20"09j |Liu et al ||20T0] >. In the figures that 
follow, we plot the fraction of isolated primaries that have 
a secondary within a given projected distance / sa t(rfproj) = 
N scc (< dproj)/A / pr i. Using the fraction rather than the raw 

6 The quantitative results here are only valid for the choice of halo and 
subhalo definitions used by MS-II (See j(3j. Other halo-finding methods may 
result in quantitatively slightly differenthalo properties, but are unlikely to 
result in different qualitative results (Kne be et al.|2011) , particularly given 
how the subhalos we consider here are isolated and well-resolved. 



number is preferred as it will be less sensitive to the effects of 
cosmic variance and the effects of cosmology on the overall 
abundance of dark matter halos. 
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FIG. 4. — Lower panel: Fraction of isolated primaries with a secondary 
within a particular projected distance for the uncorrected SDSS samples. 
The shaded regions are the uncertainties assuming Poisson statistics. The 
solid (blue) line is for the clean sample, the dashed (green) is for the any 
z, and dotted (red) is for the full sample. See Table u\ for a description of 
these samples. Upper panel: Selection function for spectroscopic sample, 
^scl = /clean//Anyz with associated Poisson errors. This selection func- 
tion is applied to the full sample for following plots. 

In the lower panel of Figure HI we plot this fraction for our 
observational samples. Solid (blue), dashed (green), and dot- 
ted (red) lines show the distributions for the Clean, Any z, 
and Full sample, respectively. See Table [T| for a description 
of these samples. We also show shaded regions around these 
lines indicating error bars assuming Poisson statistics. 

Because the SDSS spectroscopic survey has missed some 
secondaries due to failed redshifts or fiber collision, the SDSS 
spectroscopic survey suffers from incompleteness that is dis- 
joint from the biases we explored with the MS-II simulation 
in the previous section. To correct for this incompleteness we 
make a simple assumption in the SDSS data: that the true se- 
lection function relative to the full sample is the same as the 
the clean sample relative to the any z sample. We show this 
selection function F sc \ = / C iean//Anyz in the upper panel of 
Figure[4]with the associated Poisson errors; it allows us to de- 
termine the true clean fraction in the SDSS sample by taking 

/sat ~ /clean, coir = -^scl/full- 

The corrected nearest secondary fraction, /clean con •> ls plot- 
ted as a function of projected separation in the left panel of 
Figure [5] as the solid (blue) line with errors due to Poisson 
statistics shown as the blue shaded region. Over this, we plot 
the MS-II Simulated Clean sample for all 250 lines of sight as 
the thin green lines. We also show the fiber radius at the maxi- 
mum redshift of our sample as the vertical dashed (black) line. 
This radius is also the approximate physical (rather than pro- 
jected) separation at which baryonic effects not included in 
MS-II render the simulation less reliable, so we consider re- 
sults inside this distance subject to additional uncertainty that 
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FIG. 5. — Fraction of isolated primaries with a secondary within a particular projected distance in cumulative (left panel) and differential (right panel) form. 
The thick blue line on the left panel is for SDSS galaxies in the clean sample, corrected for spectroscopic selection (/doan. corr. see text), with Poisson error bars 
indic ated by the shading above and below the line. The thin green lines are the cumulative fraction of the simulated clean sample for 250 lines of sight in MS-II 
(see |3.2( . This operates as an estimate of the effect of cosmic variance. The vertical dashed (black) line indicates the distance corresponding to an SDSS fiber 
diameter at the maximum redshift of the sample, and also the approximate distance at which baryonic effects may render the MS-II sample uncertain. The right 
panel presents the same information in differential/binned form, giving the fraction per bin in d pT0 j. Points (blue) with error bars are for the SDSS clean sample 
and with Poisson errors, and the (green) bars are la ranges for the MS-II sample lines of sight. The reduced chi-squared value comparing the observations to 
simulation is also shown in this panel, and it indicates no statistically significant deviation between the two (if the bins are uncorrelated). For both panels, the 
observed fraction is corrected for selection as described in |j4]by multiplying the fraction in the full sample by F se i of Figure[4] 



is difficult to quantify. The same data are plotted in differen- 
tial form in the right panel of FigurelS] The number of SDSS 
galaxies in each d pro j bin is given by (blue) points with Pois- 
son error bars, and the MS-II galaxies are given by (green) 
bands representing la scatter over all lines of sight. From this 
plot, it is clear that the SDSS clean sample and MS-II match 
over most of the range, although with a moderately signifi- 
cant dip at ~ 80 kpc. Further, we compute the chi-squared for 
these data, and find xjLd = 0.67, with a probability p = 94% 
of a x 2 this high or higher if the simulation is a model for 
the observations. These results indicate a remarkable consis- 
tency between the SDSS observations and the predictions of 
ACDM. There is no overabundance (nor underabundance) of 
luminous satellites of isolated galaxies. 

Before moving on, we draw attention to the fact that we 
have only discussed the projected separation distribution for 
rfproj < 250 kpc. This is intentional. As mentioned in S|2j 
the targeting scheme of the SDSS spectroscopic survey intro- 
duces a bias at radii larger than this that is difficult to quan- 
tify accurately. If a background cluster (or group) is close 
enough to a primary that it is within the projected distance 
expected for a secondary, the cluster galaxies will be con- 
sidered as possible secondaries. Because there are multiple 
bright neighbors in the cluster, only one will be assigned a 
fiber by the SDSS tiling algorithm, and the others will have no 
redshift. If assigned the redshift of the primary, often at least 
one of these galaxies will fall into the magnitude range of a 
secondary. Hence, in the full sample (Table [TJ, they will be 
designated the secondary for that primary if there is no other 
nearer secondary. Thus, the background galaxy is taken to be 
a secondary, excluding a true satellite that could be present at 



a larger projected distance. This effect is not fully accounted 
for by the use of / c i ca n,corr because the number of potential 
secondaries in the no z sample increases by the number of 
cluster galaxies without redshifts, while the any z sample only 
increases by one. At large separations, this effect becomes 
significant: at d pro j > 350 kpc, many secondaries in the no z 
sample are visually associated with background cluster galax- 
ies. Fortunately, the presence of a background cluster depends 
strongly on d pro j due to the d 2 scaling of background volume. 
Visual inspection of candidate secondaries at large separation 
in the full sample shows that this effect does not become sig- 
nificant until d pro j > 250 kpc. Another reason we have cho- 
sen to limit our analysis to rf pro j < 250 kpc has to do with 



the analysis of purity and completeness described in i 3.3 the 
completeness and purity corrections conspire to be ~ 1 at that 
radius, so the fraction at 250 kpc in Figure|5]can be interpreted 
as an estimate of the true fraction of satellites within the virial 
extent of the halo. Thus, we expect that ~ 40% of isolated 
~ L.J, galaxies have a luminous (LMC-like) satellite within 
their dark matter haloQ 

4.2. Pairwise Velocity Distributions 

The line-of-sight velocities for the SDSS galaxies are di- 
rectly determined from the measured redshifts of the galax- 
ies in the SDSS sample, and they are easily measured for 
MS-II halos by projecting the halo velocity along the line-of- 
sight. In Figure 6] we show the pairwise velocity distribution 
Av = c(z pr j — z soc ) for the clean sample and the simulated 
clean sample. We compute the dispersion for the observed 

7 T he exact fraction is a function of the host halo mass, as discussed in 
i|44]andi)6] 
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FIG. 6. — Distribution of Av = c(^ pr i — z scc ) for the clean sample (solid 
blue histogram), the clean-like sample from MS-II (dashed green). The KS 
test yields pks = 33%. The pairwise velocity dispersion in the observed 
sample is a = 161 kms -1 . 

sample and find a = 161 kms -1 , of interest as a po ssible 
tracer of host or cosmological properties (e.g. |Strauss et al.| 
1998 ). The pairwise velocity dispersion also provides an indi- 
rect test of the halo mass-luminosity assignment that we have 
used in our simulation (e.g. if our halos were too massive 
then our simulated velocity distribution would presumably be 
too large). As illustrated in Figure [6] the simulated and ob- 
served distributions agree quite welL We also compute the 
Kolmogorov-Smirnov (KS) statistic between the two distribu- 
tions. This test shows that the simulation and observations are 
consistent, with a pks = 33% probability of the observed KS 
D statistic if they are drawn from the same distribution. This 
result is consistent with no velocity bias between the galaxies 
and their subhalos (e.g. Blan ton et al.|199 9), although it says 
nothing about bias of the subhalos relative to the smooth host 
dark matter distribution (e.g Colin et al. 2000 ). Thus, the rela- 
tive line-of-sight velocities in SDSS for likely bright satellites 
around L. t hosts are fully consistent with the predictions of 
ACDM. 

4.3. Luminosity Function 



In f 4. 1 we show that the MS-II mock galaxy catalogs (and 
by proxy, ACDM predictions) match the observed projected 
separation distribution of luminous satellites. Here, we com- 
pare the luminosity function implied by the observed and sim- 
ulated samples. While a luminosity function of our volume- 
limited observational sample is straightforward to compute, 
comparison to the simulation is complicated by the need to 
use abundance matching to m ap g alaxy luminosities to halos. 
We proceed as described in £3.1 but assign a luminosity to 
each halo instead of just computing mass ranges for the pri- 
mary and secondary samples. The resulting luminosity func- 
tions are shown in Figure [7] for the selection corrected clean 
sample (points) and the simulated clean sample (bars). 
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FIG. 7. — Comparison of luminosity functions. Brighter than M r < 
—20.8, the luminosity function for clean sample primaries is shown as (ma- 
genta) points, and selection corrected clean sample secondaries are shown on 
the faint end as (blue) points. Error bars on the points are Poisson. Shaded 
bars are for the 68 % bounds of the 250 realizations of the simulated clean 
sample, again for primaries on the bright end (red) and secondaries fainter 
than M r < —20.8 (green). Luminosities have been assigned to subhalos 
with abundance matching, as described in j )3.1| Also shown are reduced 
X 2 values for a \ 2 test comparing the clean and si mulated clean samples . 
Dashed black curves are the luminosity function from |Blanton et ah|^2005a) . 
The simulated clean and Blanton et al. ( 2005a| curves are scaled to match the 
normalization of the clean sample. 

We do not directly compare the normalization of the lumi- 
nosity functions, as the selection of primary galaxies for the 
SDSS spectroscopic survey results in a highly uncertain sur- 
vey volume (motivating use of the primary fraction in S4.1 1. 
Instead, we adjust the normalization of the simulated clean 
sample to match the observed clean sample and compare 
shapes. As the reduced x 2 values in Figure[7]indicate, the lu- 
minosity function for secondaries matches well within the er- 
ror bars, while the primaries are marginally deviant. This level 
of agreement is noteworthy, as subhalo abundance matching is 
not well constrained at these low masses/luminosities. While 
it is hi ghly effective where calibrated for the full galaxy sam- 
ple ( Ki-avts oveFal |2004[ |Conroy & Wechsler|2009"l |Moster| 
|et al.|2010| and references therein), abundance matching may 
not be as"er fective for biased subsets of galaxies (e.g. Tollerud 
et al.|20lT i. Thus, the disagreement for the primaries may be 
due to a preference for a particular subclass of more or less 
massive galaxies (e.g. red/passive vs. blue/starforming) to 
host secondaries. The abundance matching does not capture 
the differing luminosity functions of these two classes and 
hence would not match the obs erved sample. Alternatively, 
the bias described at the end of { 3.1 may account for the dis- 
agreement, as it would manifest as an enhancement in the ob- 
served sample in the faintest bin. Removing this faintest bin 
lowers the reduced x 2 to 1 . 1 , suggesting it may be the primary 
effect. 

The luminosity functions in the simulation and observa- 
tions agree reasonably well, further suggesting that ACDM 
matches the observations. However, this is not necessarily a 
strong test compared to that discussed in the previous sub- 
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FIG. 8. — Fraction of simulated primaries with a secondary within a given 
projected distance compared to a more general set of host halos. The dashed 
(green) line is for the simulated clean sample and its secondaries. The dot- 
dashed (black) line is the distribution of the 3D nearest object (nearest in r 
rather than d pro j) for all hosts in the MS-II SDSS-like volume. This sample 
is not subject to assembly bias effects in the way the clean sample may be, as 
it is not explicitly selected to be isolated. The thin solid lines show fractions 
for the same hosts sub-divided into quintiles in 'Umax (^max.host decreasing 
from upper-most to bottom line). The median u max host values for each line 
are indicated in the figure legend. 

section. The (blac k) dashed lines in Fig ure [7] show the lumi- 
nosity function of Blanto rTet al.| ( |2005a| l in the primary and 
secondary regimes. It is clear that the shape of these luminos- 
ity functions roughly match the clean observed and simulated 
samples. Thus, this primarily shows that the clean samples 
very roughly trace the overall luminosity function (the overall 
luminosity function implied by the simulations must match 
by construction). This is, however, a necessary consistency 
check that these samples (and hence ACDM) pass. 

4.4. Host Halo Sample Biases 

We have shown that the isolated (clean) SDSS sample 
matches the simulated clean sample from MS-II remarkably 
well. We now ask whether our isolation criteria have forced us 
to select a special set of host halos that are not representative 
of the average dark matter halo in their satellite properties. 
That is, have we measured the satellite statistics that are rep- 
resentative of all ~ L.J, halos or just a special group of isolated 
~ L.J, halos? Su ch a concern is motivated by the possibility of 
assembly b ias ( |Gao et al.|2005| [Wechsler et al. 12006} |Croton| 
|et al.|2007) . 

In Figure|8]we use the MS-II simulation to plot the fraction 
of primaries with a secondary within a particular projected 
separation for a more general sample of halos. The Simulated 
Clean sample (also shown in Figure[5]) is plotted as the dashed 
(green) line. The dash-dotted (black) line plots the separation 
distribution for all host halos in the MS-II SDSS-like volume 
in the primary mass range, regardless of whether or not they 
meet our isolation criteria. Additionally, secondaries are as- 
signed by identifying the nearest object in 3D (not projected) 
distanc e, to eliminate the projection effects described at the 
end of S3.3 We also split this sample into quintiles by host 
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FIG. 9. — Number density profile (3D) of nearest satellite. Points (blue) are 
computed as the Abel Transform of the observed selection corrected clean 
sample separation distribution, normalized to give a fraction of 42% at 250 
kpc. The error bars are generated by a combination of Poisson errors and 
bootstrap resampling (with replacement). The dashed (black) line is the best 
fit power law (p <x r~ 1,s ) for this data set, normalized to give t he same 
fraction integrated over the range shown. The dotted (green) line is a |Burkert| 
1 1995 1 profile fit to the observed data, with core radius ro = 68 kpc and 
nrj = 1.6 X 10 — 7 kpc — 3 . Finally, the solid (red) curve is the 3D density 
profile of the simulated clean sample from MS-II. 



halo u max and plot these distributions as the thin solid lines in 
Figure [8] As described in more detail in S]6] examining these 
lines clearly shows that the exact fraction at a given separation 
is strongly dependent on the choice of mass ranges used for 
the primary, and hence care must be taken comparing similar 
samples in other works to ensure they are closely-matched in 
primary /host properties. 

While there is a discrepancy of ~ 4% at 250 kpc between 
the Any "u ma x,host distribution (black dash-dotted) and the 
Simulated Clean (green dashed), this is only marginally larger 
than the scatter over multiple pointings in the mock galaxy 
catalogs, and is < ler compared to the Poisson errors in the 
observational sample. Additionally, the mass bin centered 



around the median i> n 



.host 



for the abundance-matched sam- 



ple (211 km/s; magenta line) is even closer to the Simulated 
Clean distribution. Hence, while the SDSS clean sample at 
first appears rather tuned and specific in its parameters, it is 
in fact a fairly general sample that is representative of true 
satellites of ~ L* hosts. Specifically, there is no sign that 
differences in halo assembly time associated with large-scale 
clustering have given rise to any observable bias in our iso- 
lated clean sample compared to the general halo population. 

4.5. 3D Density Distribution of Satellites 

We now turn to an estimate of the 3D distribution of satel- 
lites directly from the data themselves. We start with the 2D 
(selection-corrected, /clean, corr) clean profile from the right 
panel of Figure B] and determine the on-sky surface density 
in these annuli. From this surface density, we numerically 
perform an inverse Abel Transform to approximate the 3D 
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density profile under the assumption of spherical symmetr^ 
The resulting 3D number density profile obtained by this 
prescription is shown in Figure [9] as the (blue) points. This 
profile is by necessity coarsely binned to reduce numerical 
noise due to small number statistics, and the missing bins at 
small radii are due to the declining number of secondaries 
at radii inside the fiber radius. We also show error bars ob- 
tained from a combination of Poisson error and bootstrap re- 
sampling to simulate other choices of binning. We fit a power 
law to the data (black dotted line), normalized to integrate 
to 42% from the inner most point to 250 kpc, obtaining a 
slope of -1.8 (and normalization 6.6 x 10~ 5 ). Also shown 
is the 3D density profile of satellites of the simulated clean 
sample (green solid line), normalized to match. The den- 
sity profile in the simulations transitions from a nearly flat 
core to a profile falling off as roughly r~ 3 . This profile 
is consistent with the Abel inverted observed profile, moti- 
vating a fit to a Burkert (1995) profile as an approximation 
to the observed profile. The resulting fit has a core radius 
ro = 68 kpc, or log(ro/kpc) = 1.8, and central number den- 
sity no = 1.6 x 10~ 7 kpc~ 3 . These estimates of the profile, 
however, come with the crucial caveat that the shape (partic- 
ularly the presence of a core) cannot be truly constrained by 
the observations without a much larger data set that lacks the 
systematic biases induced by fiber collision. 

5. PHOTOMETRIC DISTRIBUTION OF SECONDARIES AND 
COMPARISON TO THE LMC 

Having compared the observational sample to the ACDM 
predictions in MS -II, we now overview the photometric char- 
acteristics of the secondaries. In this section we focus on the 
clean sample, as it is the sample that has the highest fraction 
of likely satellites. This sample is not complete due to instru- 
mental failures and the effect of distant background clusters 
(see the discussion in Sj4|. However this incompleteness is 
primarily random, and hence should not systematically bias 
the photometric properties of the sample^] 

5.1. Photometric Properties of Secondaries 

In Figure [10] we show the g — r, M r CMD for the secon- 
daries in theclean sample as the white circles, plotted over 
all galaxies within the volume limit that have magnitudes in 
the secondary range. We also show the color distribution in 
Figure[TT]for the clean sample as the shaded (blue) histogram 
and the Control sample (a collapsed version of the background 
of Figure [TO]) as the solid lined (green) histogram. The Con- 
trol sampleTwhich includes both satellites and non-satellites, 
cle arly shows the "blue c loud'T'red sequence" color bimodal- 
ity ( |Strateva et al.|[200T] l and hence below, when discussing 
blue cloud galaxies, we refer to galaxies with g — r < 0.6. 

The most striking feature of Figure 1 1 is the fact that the 
color distribution of the Clean sample secondaries is visi- 
bly different from the Control, confirmed by KS test at the 
Pks ~ 10~ 7 level. For the secondaries, 52% are on the red se- 
quence, while in the Control sample, only 38% of galaxies are 
red. This is even more surprising given that our Control sam- 
ple has galaxies in all environments, not just isolated (low- 

8 While the satellite distribution of any individual host is not necessarily 
spherically symmetric, the halos are oriented randomly relative to the ob- 
server in an isotropic universe, and thus the ensemble average for all pri- 
maries is spherically symmetric. 

9 While fiber collision may bias the sample away from very close pairs, 
these are always a small fraction of the sample, as the total volume available 
for satellites at small separations is very small. 
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FIG. 10. — Color-magnitude diagram (g — r and M r ) for the secondaries 
with d pI0 j < 250 kpc. White circles are the clean sample, and the under- 
lying hex-binned plot is a Hess-like diagram indicating the number of SDSS 
galaxies in each bin (following Figure|2j. The LMC is shown as the three star- 
shaped symbols for the t hree color esti mates: from left to right, RC3 (red), 
Eskew & Zaritsky 2011 (magenta), and |Bothun & Thom pson 1988](blue). 
The black-star shaped symbol is for M33 (see text). 

density) environments. A control sample selected to be iso- 
lated like the primaries would be significantly bluer than aver- 
age due to the density-morphology or densit y-star formation 
rate relations ( |Dressler|[T980| |G6mez et aL"1|2003| l. Because 
our Control sample is not selected in this manner, the differ- 
ence between secondaries and similar isolated field galaxies 
should be even greater than shown here. 

This result implies that star formation is suppressed even 
for single satellites of isolated host galaxies. This suggests 
that h arassment by other luminous satellites ( |Moore et al. 



1996 ) cannot be a major cause, and thus other effects such as 



strangulati on or ram press ure stripping must be at play (e.g. 
|van den Bosch et al.|2008[ and references therein). Given the 
nearly complete lack of gas for MW satelli tes other than the 
MCs, however (Grcevich & Putman 2009), the mechanism 
driving morphological transformation for these objects may 
be different from that in denser environments. A more in- 
depth analysis of the secondaries in this context is warranted 
but beyond the scope of the current work. 

5.2. Local Group Photometric Estimates 

To determine what the clean secondary sample can tell us 
about the MW/LMC system, we must estimate the integrated 
absolute magnitudes and colors of the MW and LMC as the 
SDSS would observe them from outside the LG. Our SDSS 
data set focuses on the r- and g-bands, where signal-to-noise 
is best, so below we estimate the r-band absolute magnitudes 
of the MW and LMC, as well as their g — r color. 

Determining the luminosity of the MW in the optical is no- 
toriously difficult due to the high extinction in the disk, as 
well as the necessity of precise photometry o ver a large frac- 
tion of the sky. An often-cited estimate from |van den Bergh 
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TABLE 2 

LMC Integrated Photometric Properties 



I I Clean Sample Secondaries 
■ ■ Control Sample 




1.0 



FIG. 11. — Color distribution of secondaries in the observational data set. 
The shaded (blue) histogram is for secondaries in the clean sample, while 
the line (green) histogram is for the control sample (all objects in the volume 
limit and magnitude range for secondaries). The KS test yields pks ~ 10 — 7 . 
Short dashed vertical li nes show LG galaxy color e stimates: from l eft to right, 
LMC from RC3 (red), |Eskew & Zaritsky] j20TT] magenta), and |Bothun &] 
| Thomps on (1988 blue); and M33 (black, see text). The p values in the legend 
give what percentile of the clean sample (the shaded histogram) that color 
corresponds to. The arrow for the LMC is placed at the mean of the three 
color measurements. 

(2000) is My = -20.9. To convert to SDSS r band, we deter- 
mine the average color of a sample of Sbc galaxies (the type 
for t he MW according to Ivan den B ergh 2000} in the SDSS 
from |Nair & Abraham|(|20l0|>, findin g < g - r >sbc= 0.60. 
We then use the |J ester et al.| ( |2005a| ) photometric transforms 
and convert from Vega to SDSS AB magnitudes, yielding 
M rj Mw = —21.2. This is within 0.1 mag of the median of 
the SDSS sample of the primaries (and also M*, r, by design 
of the sample). 

For the LMC, a somewhat different set of complications are 
present. Even the high surface brightness parts of the LMC 
cover a relatively large area on the sky (> 100 sq. degrees), 
and hence accurate measurements of the LMC's integrated 
magnitudes require photometry with consistent background 
subtraction over the entire area. Hence, the overall magni- 
tudes for the LMC are rather uncertain. The color is eas- 
ier to measure because it does not require consistent calibra- 
tion over the entire area. Uncertainties in the distance to the 
LMC also propagate to uncertainties in the integrated absolute 
magnitude, but the uncertainties in the surface photometry are 
much larger. Thus, we assume a distan ce to the LMC o f 50 
kpc (distance modulus 18.50, e.g. Pietrzynski et al.|2009 and 
references therein). 

With these complications in mind, we consider the fore- 
ground extinction-corrected apparent magnitude from the 



Third Reference Ca talogue of Bright Galaxies (RC3 de Vau 
|couleurs et al.|199T| l for our estimate of the LMC's r band ab 
solute magnitude. The RC3 gives U,B,V = 0.52, 0.57, 0.13, 
which we adjust to SDSS colors using the transformations of 
Jester et al.|(|2005b|), yielding M r , LM c = -18.44. The g-r 



Source 


Method 


M r 


g-r 


de Vaucouleurs et al. ( 1991 1 


Phot. Surf. Phot. 


-18.44 


0.22 


[Bothun & inompsonl 1988) 


CCD Surf. Phot. 


-18.93 


0.34 


|asKew & ZantsKyjpbu} 


Population Synthesis 


-18.44 


0.27 



color implied from the RC3 magnitudes is (g— t)lmc = 0.22. 
A second determination of the LMC integrated properties 



from surface pho tometry results from the work of Bothun & 
Thompson (1988). These observations were obtained with a 
CCD camera mounted to a standard camera lens, providing a 
wide enough field of view to cover a sizable fraction of the 
LMC in a single exposure. With foreground extincti on cor- 



rection for the LMC assuming E(B - V) = 0.075 (Schlegel 
et al. 1998 ), these observations and the above transforms yield 
(ff-f)LMC = 0.33. While the field of view was still not large 
enough to cover the luminous extent of the LMC, an extrapo- 
lation of an exponential profile yields M r lmc = —18.93. 

We also consider an alternative to surface photometry: 
stellar population synthesis. This method provides colors 
by computing the luminosity-weighted contribution of the 
present-day stellar populations with an age and metallic- 
ity distribution based on an empirical star formation history 
(SFH) and final metallicity. This method has the advantage 
of not requiring the careful calibration necessary for surface 
photometry, but has the disadvantage of bei ng more model- 
depend ent. We consider here the results of Eskew & Zaritsky 
pHl) based on the LMC SFH of |Harris & Zaritsky | (2009). 
From these stellar populations, the color of the LMC at the 
present day is (g — rj^ MC = 0-27. Adopting the stellar mass 
implied from the SFH of Harris & Zaritsky (2009) then yields 
M r LMC = —18.44. These results differ from the surface 
photometry; they do not include the effects of reddening and 
extinction internal to the LMC, as they assume all stars in the 
synthetic population contribute their full luminosity. 

Based on the above measurements, we adopt A/ r .LMC — 
-18.44 (L r = 1.75 x 1O 9 L ) from the RC3 catalog as our 
fiducial estimate for the LMC, as it is a direct measurement 
more akin to SDSS observat ions (although it is consistent 
within error bars with the Harris & Zaritsky (2009 ) result), re- 
quiring no extrapolation. While this luminosity is marginally 
fainter than the overall median of the secondaries in our SDSS 
sample (M r = —19.3), it is close to the median of the blue 



cloud galaxies apparent in the background of Figure 10 Addi- 
tionally, there is no compelling reason to favor the g — r colors 
of any of these results over the others, and thus we consider 
all of them. Thus, we summarize the photometric estimates 
for the LMC in Table 

5.3. LMC and M33 in Context 



An additional feature of Figure 10 is the location of the 
LMC using our three color estimates and the RC3 magnitude 
as the star-shaped symbols. We also indicate the LMC color 
estimates as the left three vertical lines in Figure [TT] For all 
colors we consider here, the LMC is remarkablyBlue rela- 
tive to other satellites in its ma gnitude range, in the 0.3 , 1 .3 , 
and 9.9 percentile for the RC3, Eskew & Zaritsky (201 1}, and 
Bothun & Thompson ( 1988 ) colors, respectively (also shown 
in the legend of Figure 11). In the particular case of the RC3 
color, the LMC is the bluest object in the sample. While the 
other colors are somewhat less extreme, they still place the 
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LMC firmly on the blue side of the blue cloud. For the more 
luminous satellites, this placement is partly due to the fact that 
the red sequence becomes more significant at brighter magni- 
tudes, but even restricted to the magnitude range close to the 
LMC, the LMC is unusually blue. Thus, while the LMC/MW 
system is not particularly unusual in its existence from the 
perspective of the sizes of the galaxies involved, it is very un- 
usual in the sense that the LMC has an anomalously high star 
formation rate. 

This remarkably large star formation rate for the LMC, es- 
pecially when compared to similar satellite systems, begs the 
question of its cause. The most straightforward explanation 
for this effect is that the LMC is currently undergoing trig- 
gered star form ation due to tidal or gasdynami cal perturbation 
by the MW ([Kennicutt etal.|1987 ; Funes et al. 20041 |Woods 



[etaT1|2006l |Woods & G eller 2007). This is a 



so consistent 



with starformation histories o f the LMC (Harris & Zaritsky 
|2009| [Maschberger & Kroupa|2010| l that indicate an increase 
in the star formation rate in the last 100 Myr - 1 Gyr (although 
there are some earlier possible peaks). Thus, the very blue 
colors of the LMC relative t o similar objects are consistent 
with recent proper motions (Kallivayalil et al. 2006; Piatek 



et al.||2008) that place the LMC on first infall JBesla et al. 



2007| |Boylan-Kolch in et al.|20l0 a). Together with the result 
that other luminous satellites are red and all MW satellites 
aside from the MCs are no n-starforming and extremely gas- 
poor ( |Grcevich & Pu tman 2009 ), the complete picture sug- 
gests that this star formation rate cannot be sustained for very 
long. 

Finally, we consider the LG galaxy M33. If placed in the 
SDSS, M33 would appear as a secondary of M31 p] with a 
projected sep aration of 208 kpc (ass uming an M31 distance of 
791 kpc from |Fiorentino et al.|2010| > and Av < 500 kms~\ 
If we take M33 to have My = -18.87 and B - V = 0.65 
following van den Bergh (2000), converting to SDSS bands 
as described above for the MW gives M ri M33 = —19-06 and 
g — r = 0.44. Thus, in a cosmological context, Figure [3] 
gives an ~ 86% probability that M33 is actually a satellite 
of M31. However, its g — r color is quite typical for star- 
forming secondaries (see Figures 10 and [TT) , suggesting it 
does not have the significant triggered star formation evident 
in the LMC. 

6. COMPARISON WITH OTHER WORK 

We now consider these results in the context of r ecent work 
examining the frequency of MC-like satellites. Liu et al. 
(2010) examined a similar set of MW/LMC system analogs 
in the SDSS photometric sample and determined the abun- 
dance of satellites. Our selection is distinguished from these 
by the use of spectroscopic redshifts to define our clean sam- 
ple of secondaries, somewhat different selection criteria, and 
the use of separation distributions instead of a single distance 
(they adopted 150 kpc). While our use of the spectroscopic 
sampl e reduces our vo lume relative to that of the main sam- 
ple of |Liu et aL] ( |2010[ ), individually, our objects have a much 
higher likelihood of being true satellites because our redshifts 
are spectroscopic rather than photometric. Nonet heless, our 
results are consistent within the error bars - the ILiu et aLl 



(2010) result of roughly 20% probability to host one satellite 
matches our result at 150 kpc. We note that, as described in 

10 M31 would be considered isolated from a distant viewer's perspective 
unless the MW-M31 separation vector is close to the line of sight, which 
would occur ~ 35% of the time for random orientations. 



j |3.3| 250 kpc is more appropriate for asking the specific ques- 
tion of whether or not a host has a satellite within its virial 
radius. 

IBoylan-Kolchin et al.| ( [2010b| ) and |Busha et~aT1d2010a| ) ex- 
amined exactly this question using high resolution n-body 
simulations (Millenium II and Bolshoi, respectively). They 
both found similar results, finding ~ 5 — 30% of MW analogs 
host LMCs. While this may appear at first glance to be at dis- 
cord with our results, examination of Figure [8] explains this 
apparent discrepancy. As the thin solid lines show, the frac- 
tion of hosts with satellites within a given separation depends 
quite stron gly on the halo mas s. This fact is taken ad vantage 
of by |Busha et al.| ( [20T0a| and |Boylan-Kolchin et aT] ( |2010a i 
to constrain the properties of the MW based on the existence 
and observed prope rties of the MCs. Further, as illustrated in 
Busha et al. (2010a), the satellite mass range also has a (some- 
what weaker) effect on the satellite fraction. Thus, while the 
exact definition of both "MW-like" and "LMC-like" will af- 
fect the exact fraction, the effect is generally only a factor of a 
few for reasonable definitions. Further, when we test match- 
ing definitions, we find identical answers within the relevant 
error bars. 

Additionally, in the context of the color distribution for the 
SDSS clean sample relative to the LMC, the low num bers of 
LMC-like satel lites reported in |James et al. (2008 ) and James 
& Ivory (2010) become clear. Both are ifa-selected surveys 
searching MW-like hosts for satellites with LMC-like R-band 
magnitudes and Ha emission (indicating significant star for- 
mation). Both report LMC-like satellites are rare, implying 
the LG is a statistical outlier. Figure 1 1 implies that this ap- 



parent anomaly is in fact due to the over-abundance of r ed sec- 
ondaries relativ e to blue for the primaries in question. |James 



& Ivory (2010 1 Figure 5, in particular, indicates that 8.2% 
of hosts have star-forming satellites in the "LMC" category 
within 100 kpc (including their specified volume correction 
of 1.96). While the exact selection criteria are slightly dif- 
ferent, this result is fully consistent within Poisson error bars 
of the fraction of hosts with blue cloud galaxies that we find 
within 100 kpc. Further, their result that the i/a-inferred star 
formation rate of the LMC is higher than most of the other 
star-forming satellites in their survey lends further credence 
to the claim that the LMC has an unusually large star forma- 
tion rate for its status as a satellite of an isolated ~ host. 

7. CONCLUSIONS 

In this study, we examine a spectroscopic sample of galaxy 
pairs from the SDSS selected to be composed of an isolated 
~ L.J, (M r < —20.8, median -21.3) galaxy with a ~ 0.1 L* 
(-18.3 > M r > -20.8, median M r = -19.3) nearest satel- 
lite, loosely resembling the MW/LMC system. We do this by 
identifying objects with no ~ galaxies nearby (on-sky), 
and identify fainter galaxies within 250 kpc that have consis- 
tent redshifts as likely satellites. We also examine identically 
selected halos and subhalos in a mock galaxy catalog from the 
MS-II cosmological simulation. The high purity of the spec- 
troscopic sample allows us to directly compare secondaries 
to subhalos, confirming our interpretation of the observations 
in a ACDM context. It also allows for direct examination of 
the secondaries' photometric properties. This is therefore a 
cosmologically significant sample of "small-scale structure". 
From this work we draw the following conclusions: 

1. The projected radial distribution (Figure[5]l and line-of- 
sight velocity distribution (Figure [6} orbright secon- 
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daries matches ACDM predictions. This implies that 
there is neither an overabundance nor an underabun- 
dance problem for luminous (M r < —18.3) satellites 
of isolated ~ (M r ~ —21.3) galaxies. Specifi- 
cally, with minimal assumptions (abundance matching) 
ACDM can reproduce the observed properties of the 
Universe for galaxies as faint as L ~ 1O 9 L on veloc- 



ity scales of > 
kpc. 



95 kms and spatial scales of ~ 50 



2. Bright secondaries are significantly redder than a gen- 
eral sample of galaxies of the same luminosity (Figures 
[T0|and[TT|. More than half of the ~ 0.1 L* secondaries 
inhabit the red sequence, compared with only ~ 30% 
in our control sample. This suggests that environmen- 
tal quenching is operating within the halos of ~ 
galaxies, in the absence of harassment from compara- 
bly bright satellite counterparts. 

3. MW-size galaxies hosting LMC-size satellites are rea- 
sonably common in the SDSS sample. We estimate 
that ~ 42% of ~ galaxies host a ~ 0.1 satel- 
lite within the virial radius of their dark matter halos. 
Further, 10 % of isolated ~ galaxies have a bright 
secondary within 50 kpc, the distance to the LMC. The 
exact fraction depends on the mass range used for the 
host, however, fully accounting for the mild difference 
between this and other similar measurements. In any 
case, this implies that the presence of a satellite with 
the luminosity of the LMC in a MW-like galaxy dark 
matter halo is not a major anomaly in a cosmological 
context. 

4. The LMC is remarkable insofar as it is one of the bluest 
satellites in the sample. This color is a consequence of 
an unusually high star formation rate, and studies of the 
LMC must be cautious when generalizing its star for- 
mation properties to other galaxies. This also suggests 
the LMC may be undergoing strong triggered star for- 
mation if it is on first i nfall, consistent with the most 

recent pro per motions (Kallivaya lil et al.| 2006; Piatek 

et al.|2008) . 

In the end, the observations presented here represent an- 
other test passed by ACDM, and place the Local Group firmly 
in the cosmological context of the local universe. They do, 
however, beg some questions. First, at what scale, if any, do 
the standard ACDM abundance matching assumptions break 
down? This work shows that it is successful for u max > 



95 km s , while the results of Boylan-Kolchin et al. ( 2011 
suggest problems with MW satellites for v max ~ 50 km s~ . 
This implies there may be a break down scale somewhere in 
the range 50 < v max < 95 kms -1 , but there are no clear 
causes of such a break down in this range. An additional 
important question is that of what causes the quenching of 
star formation in these secondaries. Finally, is the remark- 
able blueness of the Milky Way's LMC just chance timing, or 
does it imply something more fundamental? The answers to 
these questions are well beyond the scope of this work, but it 
is clear that these isolated host/satellite pairs, while directly 
comparable to the wider ACDM cosmology, may teach us as 
much about our own small corner of the Universe. 
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